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ABSTRACT 

We have used the latest HI observations of the Small Magellanic Cloud (SMC), obtained with the 
Australia Telescope Compact Array and the Parkes telescope, to re-examine the kinematics of this dwarf, 
irregular galaxy. A large velocity gradient is found in the HI velocity field with a significant symmetry in 
iso- velocity contours, suggestive of a differential rotation. A comparison of HI data with the predictions 
from tidal models for the SMC evolution suggests that the central region of the SMC corresponds to 
the central, disk- or bar-like, component left from the rotationally supported SMC disk prior to its last 
two encounters with the Large Magellanic Cloud. In this scenario, the velocity gradient is expected as 
a left-over from the original, pre-encounter, angular momentum. We have derived the HI rotation curve 
and the mass model for the SMC. This rotation curve rapidly rises to about 60 km s~^ up to the turnover 
radius of 3 kpc. A stellar mass-to-light ratio of about unity is required to match the observed rotation 
curve, suggesting that a dark matter halo is not needed to explain the dynamics of the SMC. A set of 
derived kinematic parameters agrees well with the assumptions used in tidal theoretical models that led 
to a good reproduction of observational properties of the Magellanic System. The dynamical mass of 
the SMC, derived from the rotation curve, is 2.4 x 10^ M0. 
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1. INTRODUCTION 

The Small Magellanic Cloud (SMC) is a nearby^, gas- 
rich, dwarf irregular galaxy. Its morphology, dynamics, 
and evolution are highly complex, and must have been 
heavily influenced by gravitational interactions with the 
nearby Large Magellanic Cloud (LMC) and the Galaxy 
(Putman et al. 1998). As a dwarf irregular galaxy, the 
SMC is different from our own Galaxy in many respects, 
having a lower heavy element abundance, significantly 
lower dust content, and a consequently stronger interstel- 
lar radiation field (Stanimirovic et al. 1999; Stanimirovic 
& Lazarian 2001). 

Many recent high resolution studies have shown that 
dwarf galaxies in general have a very dynamic interstellar 
medium (ISM), structured mainly by star formation and 
its aftermath. The most obvious example of star- formation 
creativity are numerous expanding shells of gas (Puche 
et al. 1992; Staveley-Smith et al. 1997; Kim et al. 1999; 
Stanimirovic et al. 1999; Walter 1999). Dwarf galaxies 
are also usually found to have large, and dynamically im- 
portant, halos of dark matter (Mateo 1998). The relative 
roles of a galactic rotation, pressure support, dark matter, 
and magnetic fields on the 3-D structure and dynamics 
of these galaxies are still open questions. Furthermore, it 
is not clear yet whether different stages of galactic evo- 
lution are marked by fundamentally different phenomena, 
and whether these phenomena are selective with respect 



to spatial scales. 

It has been known for some time that the morphol- 
ogy and kinematics of the SMC traced by different stel- 
lar populations show very different properties (Gardiner 
& Hatzidimitriou 1992). This has been one of the major 
reasons to start thinking about forces other than gravita- 
tional that may have played a significant role in the forma- 
tion and evolution of the whole Magellanic System. The 
questions of morphology and kinematics are closely related 
to the long-standing and greatly controversial questions of 
the SMC's 3-D structure and depth along the hne-of-sight. 
As the results of several new optical and near-IR surveys 
are becoming available (Zaritsky et al. 2000; Cioni et al. 
2000; Maragoudaki et al. 2001), revealing structural evolu- 
tion of the SMC, as well as new N-body simulations of the 
dynamical evolution of the SMC (Yoshizawa & Noguchi 
2003), we find that it is timely to re-examine the morphol- 
ogy and kinematics of the SMC as traced by the neutral 
hydrogen (HI). 

The mapping of HI in the SMC has a long and pro- 
ductive history. After the pioneering work by Kerr et al. 
(1954) and Hindman et al. (1963), Hindman (1967) was 
the first to notice the velocity gradient in the SMC and 
to model its rotation curve, followed by Bajaja & Loiseau 
(1982). However, the velocity field of the SMC is far from a 
simple text-book example. The HI profiles, often complex 
and with multiple peaks, have caused much controversy 
in the past, being interpreted as due to either expanding 



^ Throughout this paper we assume a distance to the SMC of 60 kpc (Westerlund 1997). 
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shells of gas or spatially separate systems (Hindman 1967; 
Mathewson et al. 1988; Martin et al. 1989). An improve- 
ment by a factor of ten in spatial resolution, relative to 
previous surveys with the Parkes telescope, was achieved 
in the HI survey by Staveley-Smith et al. (1997), using the 
Australia Telescope Compact Array (ATCA). These new 
high resolution data showed that much of the HI profile 
complexity lies in the huge number (~ 500) of expand- 
ing shells. The new high resolution HI observations were 
complemented with new low resolution observations ob- 
tained with the Parkes telescope to provide information 
over a continuous range of spatial scales from 30 pc to 4 
kpc (Stanimirovic et al. 1999). 

The aim of this paper is to re-examine the HI kinemat- 
ics of the SMC, as viewed from high resolution observa- 
tions, and compare it with predictions from tidal theoreti- 
cal models. We start by summarizing the HI observations 
in Section 2. The morphology and kinematics of the SMC 
from the HI distribution, as well as viewed using other 
tracers, are discussed in Section 3. The SMC 3-D structure 
and line-of-sight depth are reviewed briefly in Section 4. 
We then compare HI data with predictions from several 
tidal models in Section 5. The rotation curve and mass 
model of the SMC are investigated in Section 6. Finally, a 
summary and concluding remarks are given in Section 7. 

2. HI DATA 

The HI in the SMC was observed with the ATCA^, a 
radio interferometer, in a mosaicing mode (Staveley-Smith 
et al. 1997). Observations of the same area were obtained 
also with the 64-m Parkes telescope. The two sets of ob- 
servations were then combined (Stanimirovic et al. 1999), 
resulting in the final HI data cube with angular resolu- 
tion of 98 arcsec, velocity resolution of 1.65 km s~^, and 
1-fT brightness-temperature sensitivity of 1.3 K to the full 
range of spatial scales between 30 pc and 4 kpc. The area 
covered with these observations is RA 00^30™ to 01^30™ 
and Dec -71° to -75° (J2000), over a velocity range of 90 
to 215 km s~^. For details about the ATCA and Parkes ob- 
servations, data processing, and data combination (short- 
spacings correction) see Staveley-Smith et al. (1997) and 
Stanimirovic et al. (1999). 

3. STRUCTURE AND KINEMATICS OF THE SMC 

3.1. HI distribution 

Fig. 1 shows the integrated HI column-density distribu- 
tion. The large-scale HI morphology of the SMC is quite 
irregular and does not show obvious symmetry. The most 
prominent features are the elongation from the north-east 
to the south-west and the V-shaped concentration at the 
east (RA 01^ 14"^, Dec -73° 15', J2000). These are usu- 
ally referred to as the 'bar' and the Eastern Wing, al- 
though their dynamical importance is still not well un- 
derstood. A 'bridge' of gas appears to connect the 'bar' 
and the Wing (note that this is not the Magellanic Bridge 
that connects the SMC and the LMC), while the arm- like 
extension of the 'bar' towards the north-east is also promi- 
nent. On smaller scales and looking at different velocity 
channels (see Stanimirovic et al. 1999), the HI distribu- 
tion appears very complex and frothy, being dominated 

^ The Australia Telescope is funded by the Commonwealth of Australia 



by numerous expanding shells, filaments and arcs. The 
total estimated mass of the HI, after correction for self- 
absorption, is 4.2 X 10* Mq (Stanimirovic et al. 1999). 
A summary of radio and optical properties of the SMC is 
given in Table 1. 

The HI column density distribution is super-imposed on 
the V-band optical image of the SMC, kindly provided to 
us by M. Bessell, in Fig. 2. In general, in the 'bar' and 
the Eastern Wing the stellar and HI distributions correlate 
well, however the HI is significantly more extended towards 
the south-west, the south-east, and the north-west. 

3.2. HI kinematics 

The intensity-weighted mean velocity along each Une- 
of-sight, or the first moment map, was derived from the 
full resolution HI data cube (corrected for short spacings) , 
and is shown in Fig. 3. A large velocity gradient, from 
91 km in the south-west to 200 km s~^ in the north- 
east, can be seen. The iso-velocity contours show some 
large-scale symmetry suggestive of differential rotation for 
the main gaseous body of the SMC. However, clear distor- 
tions are visible in the north-west, most likely correspond- 
ing to several shells and filamentary features aligned into a 
chimney-like structure at velocity ^ 123 km s~^ (Fig. 2 in 
Stanimirovic et al. 1999), and in the south-east, towards 
the Eastern Wing region, where again a supergiant shell 
was found (494A; see Stanimirovic et al. 1999). These per- 
turbations form a S-shaped feature perpendicular to the 
direction of the main velocity gradient (as shown in Fig. 3 
with a dashed line). 

The second moment map or the intensity-weighted ve- 
locity dispersion is shown in Fig. 4. This velocity disper- 
sion varies from ~5 to 40 km s~^ across the SMC. The re- 
gions with higher dispersion appear to be associated with 
the positions of the three largest supergiant shells (SGSs). 
The region around RA Ol"^ OO'^, Dec -71°30' (J2000) has 
the lowest dispersion with a mean value of ^10 km s"-*^. 
As this region corresponds to the receding hemisphere of 
the SGS 304A (Stanimirovic et al. 1999), the lower veloc- 
ity dispersion may be explained with the fact that most 
of the approaching hemisphere of SGS 304A is missing at 
the northside. 

We have investigated two effects that could significantly 
influence the observed velocity field shown in Fig. 3: (i) 
the effect of the proper motion of the SMC, and (ii) the 
effect of multiple-peak velocity profiles, 
(i) The binary motion of the SMC around the Galaxy, 
combined with its large angular extent, could have a sig- 
nificant contribution to the observed velocity field, as was 
shown for the case of the LMC by Luks & Rohlfs (1992) 
and Kim et al. (1998). This contribution consists of the 
projection of the transverse velocity of the SMC's center 
of gravity on the line-of-sight at each position. In the case 
of a large angular extent on the sky, this projection can 
vary greatly across the field of observation and cause a 
significant change in the observed velocity field. To cor- 
rect for this effect, we used values for the SMC's proper 
motion and the heliocentric transversal velocity predicted 
by numerical simulations by Gardiner et al. (1994), which 
are in agreement with the estimates based on a combi- 
nation of proper motion measurements by Kroupa et al. 

for operation as a National Facility managed by CSIRO. 
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Table 1 

Radio and optical properties of the SMC. 



Property Value 

Radio: 

Measured: 

RA (J2000)° 01^ 05™ 

Dec (J2000)° -72° 25' 

Systemic velocity'', Vgys (Gal.) 24 km s~^ 

Systemic velocity*', Vsys (Hel.) 160 km s~^ 

HI mass^ Mhi 4.2 x 10^ Mg 
Model- dependent: 

Inclination'^, i 40°± 20° 

Kinematic Position angle'', PA 40° 

VLx 60 km s-i 

Rt.^ 2.5 - 3 kpc 

M^y„ 2.4 X 109 Mg 

Optical: 

Extinction'^, i^B-v ~ 0.05-0.25 mag 

Total visible magnitude-^, Vt +2.4 mag 

Absolute visible magnitude-^. My —16.5 

Visible luminosity-^ , Ly 3.1 x 10^ Lg 

Mh/Lv 1-4 Mg/Lg 

{B - V)% 0.41 

[U - B)^ -0.23 



"The apparent kinematic center derived as the intersection of the principal kinematic axes. 
''Systemic velocity of the apparent kinematic center derived from the velocity field. 
■^Stanimirovic ct al. (1999). 
''Derived from the tilted ring analysis in Section 6. 
■^Zaritsky ct al. (2002). 
■'^de Vaucouleurs et al. (1991). 



Fig. 1. — [fl.jpeg] An HI column-density image of the SMC. The grcy-scalc intensity range is to 7 X lO^"' atoms cm ^ with a linear 
transfer function. The maximum HI column density, 1.43 X 10^2 atoms cm-2, is at position RA 00''47™33=, Dec -73°05'26" (J2000). 



Fig. 2. — [f2.jpeg] The HI column density image super-imposed on the V-band CCD image. The contour levels range from 2.0 x 10^" to 
1.3 X 10^^ atoms cm~^, with a contour interval of 3.8 x 10^" atoms cm~^. The V-band image was obtained by Mike Bessell using the 16" 
telescope at the Mount Stromlo Observatory. 



Fig. 3. — [fS.jpeg] The intensity-weighted mean heliocentric velocity field of the SMC, derived from the combination of the ATCA and 
Parkcs data. The grey scale range is 120 to 200 km s~^ with a linear transfer function. The contour levels range from 120 to 200 km s^^, with 
a contour interval of 10 km s~^. All pixels with a corresponding HI column density lower than 9 X 10"^" atoms cm"-^ were masked out. The 
direction of the largest velocity gradient is shown with a solid line, and the direction of the largest velocity distortions, forming a S-shaped 
features, is shown with a dashed line. 
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RA (J2000) 

Fig. 4. — The velocity dispersion of the SMC from the second moment analysis obtained from the combined ATCA and Parkes data cube. 
The grey-scale range is 6 to 37 km s~^. Positions of the three largest shells (37A, 304A and 494A) from Stanimirovic et al. (1999) are overlaid. 
All pixels with a corresponding HI column density lower than 7 x 10^" atoms cm~^ were masked out. 



(1994) and Kroupa & Bastian (1997). At any position 
across the SMC the proper-motion corrected heliocentric 
velocity was obtained by subtracting the projection of the 
heliocentric transverse velocity of the center of the SMC 
on the linc-of-siglit from the observed radial heliocentric 
velocity. The values for the conversion between the ob- 
served and the proper-motion corrected heliocentric veloc- 
ity range between —23 km s~^, in the north-west, and 28 
km s~^, in the south-cast. The correction for the proper 
motion, however, did not make as big difference as in the 
case of the LMC (Luks & Rohlfs 1992). The reason is that 
the angular extent of the SMC is less than half that of 
the LMC. In addition, the direction of the SMC's motion 
is perpendicular to the observed velocity gradient result- 
ing in the gradient being almost unaffected by the proper 
motion. 

(ii) Velocity profiles in the SMC are usually very com- 
plex, having double or even multiple-peak components. 
Since the intensity-weighted mean velocity (Fig. 3) is bi- 
ased towards the velocity component with higher inten- 
sity, it is not necessarily the best representation of com- 
plex velocity profiles. To avoid this bias, we found, for 
each linc-of-sight, the minimum and maximum velocity 
for which the intensity is 5% of the peak value for that 
line-of-sight. The mean of the minimum and maximum 
velocity, determined in this way, was then taken to be our 
new velocity estimate. Since such an analysis requires a 
high signal-to-noisc ratio, we only used the low resolution 
Parkes data cube. The resultant velocity field has a more 
regular 'spider' pattern than the mean one, however ir- 
regularities are still present in the north and close to the 
Eastern Wing region. We also investigated the existence 
of a bimodal velocity field, by modeling each velocity pro- 
file in the Parkes data cube with one or a superposition 
of two independent Gaussian functions. As a result, two 
separate velocity structures were depicted, with central ve- 
locities corresponding to the low and high velocity entities 



of 137 km s^^ and 174 km s^^, respectively. The position- 
velocity diagrams show almost parallel velocity fields of 
these two components. However, in many cases the ve- 
locity fields intersect stiggcsting that two separate velocity 
components may be a consequence of the statistical data 
handling of complex HI profiles with no substantial phys- 
ical meaning. 

The proper-motion corrected velocity field in the Galac- 

toccntric velocity reference frame, derived from the Parkes 
HI data cube using the 5% cut-off described in (ii), is 
shown in Fig. 5. A clear velocity gradient is noticeable 
from —15 km s~^, in the south-west, to 53 km s~^, in the 
north-east. The positions of the three largest shells in the 
SMC arc also overlaid in Fig. 5, suggesting that super- 
shells 304A and 494A may still be responsible for some of 
the perturbations in the velocity field. Some perturbations 
are also visible in the north-west. 

3.3. Comparison with other tracers 

3.3.1. SMC morphology from stellar populations 

It has been known for some time that young and old 
stellar populations in the SMC have significantly differ- 
ent spatial distributions (e.g. Gardiner & Hatzidimitriou 
(1992)). Several recent stellar surveys have especially en- 
hanced our knowledge about the different morphological 
properties of different stellar populations. The recent stel- 
lar survey of the c;cntral SMC (4°x 4°) by Zaritsky et al. 
(2000) showed that while young stars (with ages ^ 200 
Myr) exhibit an irregular distribution, similar to that seen 
in HI, the old stellar population (with ages ^ 1 Gyr) shows 
regular, undisturbed, and almost spheroidal distribution. 
Similar results were reached from the DENIS survey by 
Cioni et al. (2000) based on IJKs IR bands. These stud- 
ies also pointed out that the 'bar' and the Eastern Wing 
region may not be distinct dynamical entities, but could 
be a result of recent hydrodynamic interactions between 
the LMC and SMC's gaseous components. Maragoudaki 
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RA (J2000) 

Fig. 5. — The galactoccntric velocity field of the SMC, from tfic Parltes HI observations, witfi angular resolution of ^ 15 arcmin, corrected 
for the proper motion. The grey-scale range is —15 to 50 km s~^with a linear transfer function. The contour levels range from —20 to 50 
km s~^, with an interval of 5 km s~^. The three largest supershells (37A, 304A and 494A) from Stanimirovic et al. (1999) are superimposed. 



et al. (2001) further investigated the dynamical origin of 
the 'bar' by using isodensity contour maps of stars with 
different ages. They found similar results concerning old 
stellar populations. In addition, their data show that the 
'bar' and the Eastern Wing region were already promi- 
nent features some 0.3-0.4 Gyr ago; however, it is not 
clear whether this is a temporary extended region of star 
formation or a genuine, dynamical feature. 

3.3.2. SMC kinematics from stellar populations 

The intermediate age and old stellar populations, traced 
by planetary nebulae (Dopita et al. 1985) and carbon stars 
(Hardy et al. 1989; Hatzidimitriou et al. 1997; Hatzidim- 
itriou 1999; Kunkel et al. 2000), do not show signs of rota- 
tion. The red horizontal branch clump stars in the north- 
cast show a velocity gradient of 70 km which was in- 
terpreted as due to a distance spread of 10 kpc along the 
line-of-sight (Hazidimitriou et al. 1993). The young stel- 
lar population, traced by HI shells, shows a clear velocity 
gradient from the south-west to the north-east (Staveley- 
Smith et al. 1997), while Cephcids may have a similar 
trend with velocity (Mathewson et al. 1988). All popula- 
tions show similar mean heliocentric velocity and velocity 
dispersion regardless of age and location (Hatzidimitriou 
1999; Staveley-Smith et al. 1997). 

From an analysis of 150 carbon stars, Kunkel et al. 
(2000) estimated the inclination of the SMC orbital plane 
toi = 73° ±4° relative to the plane of the sky. They found 
that carbon stars are associated with more negative veloc- 
ity and denser HI, while low density HI seems to be de- 
void of carbon stars. This suggests that non-gravitational 
forces must be acting on the gaseous component to pro- 
vide its separation from stellar systems. Similar estimates 
for the SMC's inclination were reached using Cepheids by 
Caldwell & Coulson (1986), i = 70° ± 3°, and by Groe- 
newegen (2000), i = 68° ± 2°. 

4. THE 3-D STRUCTURE OF THE SMC 



4.1. Previous observational models 

The 3-D structure of the SMC has been a matter of 
great controversy in the past. Since early mapping of HI, 
complex HI profiles pointed to the existence of unusual 
motions of the gas (Hindman 1967). Hindman (1967) also 
suggested the first model of the SMC, as a flattened disk 
with three supergiant shells in the main body. 

From the analysis of the radial velocity distributions of 
HI, stars, HII regions, and planetary nebulae, Mathew- 
son & Ford (1984) and Mathewson, Ford & Visvanathan 
(1986, 1988) suggested the 'two separate entities' model 
for the SMC, whereby the SMC is, along most of its angu- 
lar extent, broken into two velocity subsystems (the Small 
Magellanic Cloud Remnant. SMCR, and the Mini Magel- 
lanic Cloud, MMC). In this model, the SMCR and MMC 
are separated in velocity by ~ 40 km s~^, and have their 
own nebular and stellar populations. The reason proposed 
for this disruption of the SMC was its close encounter with 
the LMC some 2x10*^ yr ago. A slight modification of the 
'two separate entities' model was suggested by Torres & 
Carranza (1987), and Martin et al. (1989), who found four, 
instead of two, different velocity components. 

Caldwell & Coulson (1986) also used Cepheids to sug- 
gest the 'bar and three arms' geometrical model of the 
SMC. The central bar is very elongated (5-to-l) and is 
seen edge-on. A distant arm of material is pulled from the 
center of the SMC to the west, while two sides of the bar 
located in the north-east and the south-west are consid- 
ered as two near arms. 

4.2. Depth along the line-of-sight 

The question of the 3-D structure of the SMC is closely 

related to the long-standing and greatly controversial issue 
of the SMC's depth along the line-of-sight. By measur- 
ing distances and radial velocities of 161 Cepheids, Math- 
ewson et al. (1986) found a great depth along the line- 
of-sight of -^30 kpc, and a distance gradient from the 
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north-east to the south- west. The large depth of the 
SMC was supported later by Mathewson et al. (1988), 
using Cepheids again, and by Hatzidimitriou & Hawkins 
(1989) and Hatzidimitriou et al. (1993) from a study of 
the intermcdiatc-agc (halo) population. 

Welch et al. (1987) have pointed out that the determina- 
tion of the Cepheid distances in Mathewson et al. (1986) 
suffers from a few possible problems: an inconsistent cor- 
rection to mean magnitude, the assumption of zero intrin- 
sic scatter for the period-luminosity (P-L) relation, and 
sample inhomogeneity. Welch et al. (1987) concluded that 
the SMC does not extend in depth beyond its tidal radius 
(4-9 kpc). Martin et al. (1989) confirmed that the young 
stars in the SMC lie within a depth of < 10 kpc. On the 
other hand, Groenewegen (2000) used near-infrared obser- 
vations of Cepheids, which are less affected by reddening 
and metallicity, and derived the depth of the SMC of 14 
kpc, assuming that the P-L relation has the same scatter 
for both the LMC and the SMC. 

Although Cepheids are used as the main distance indica- 
tors, the P-L relation requires an independent zero-point 
calibration, which is very difficult to achieve, and could 
suffer from a significant metallicity dependence. Another 
source of error in the distance determination, that has not 
been previously fully appreciated, is differential extinction 
toward the SMC. While reddening of Cepheids in the SMC 
has been extensively studied (see Welch et al. (1987) for 
a summary), a constant value has been often applied for 
distance determination. For example, Mathewson et al. 
(1986) found no sign of differential reddening and assumed 
a uniform correction of 0.06 mag for their sample of 161 
Cepheids. 

Recently, Zaritsky et al. (2002) determined the extinc- 
tion map across the SMC, using the Magellanic Clouds 
Photometric Survey, and showed that extinction varies 
both spatially across the SMC and with stellar population. 
They found that young, hot stars have an average extinc- 
tion 0.3 mag higher than old, cool stars, and that young 
stars show a significant increase in extinction along the 
main SMC ridge, from the north-east to the south-west. 
Zaritsky et al. (2002) foimd E^-v ^ 0.05 - 0.25 mag. For 
comparison, the mean value for the interstellar reddening 
previously measured by Caldwell & Coulson (1985) was 
Eb-v = (0.054 ± 0.021) mag, while Sasselov et al. (1997) 
found a sHghtly higher value of Eb-v = (0.125 ± 0.009) 
mag. 

The differential extinction across the SMC has a sig- 
nificant influence on the distance estimate. Differential 
reddening toward Cepheids in the south-west part of the 
SMC is ~ 0.2 mag higher than toward in the north-east. 
This results in the interstellar infrared absorption being 
larger in the south-west by ~0.4 mag. Hence, the Cepheid 
distances estimated by Mathewson et al. (1988) could be 
over-estimated by up to 18%. This corresponds to about 
10 kpc at a distance of 60 kpc, and suggests that the cor- 
rection for the interstellar absorption can significantly in- 
fluence the distance determination, and easily bring the 
depth of the SMC within its tidal radius (4-9 kpc). 

5. COMPARISON WITH THEORETICAL MODELS 

There are two families of theoretical models, based on 
either a tidal or a ram pressure scenario, that try to re- 



produce observational features in the Magellanic System, 
caused by interactions between the SMC, LMC, and the 
Galaxy. The 3-D structure and kinematics of the SMC 
were particularly addressed in tidal models by Gardiner 
et al. (1994), Gardiner & Noguchi (1996), and Yoshizawa & 
Noguchi (2003). The recent work by Yoshizawa & Noguchi 
(2003) provides, to date, the most detailed model for the 
evolution of the SMC. We summarize here results from 
these tidal simulations and their major predictions con- 
cerning the structure and kinematics of the SMC. 

It is important to note that none of the theoretical 
models so far predicts the bimodal velocity distribution 
throughout the main gaseous body of the SMC. This gives 
support to the idea that most, if not all, of the observed 
line-splitting comes from the combined effects of numerous 
expanding shells (Stavclcy-Smith et al. 1997). 

5.1. Model predictions 

The above three models come from N-body simula- 
tions of the gravitational interactions in the Galaxy-LMC- 

SMC system to reproduce the observed gas distribution 
in the Magellanic System, primarily those of the Magel- 
lanic Stream and Bridge. While Gardiner et al. (1994) 
modeled the SMC as a single-component disk-like system, 
in simulations by Gardiner & Noguchi (1996) the SMC 
was represented by a two-component particle system con- 
sisting of a nearly spherical halo and a rotationally sup- 
ported disk, with a disk-to-halo mass ratio of 1:1. This 
high mass ratio resulted in the original disk quickly becom- 
ing unstable, and being transformed into a bar-like struc- 
ture. The tidal model by Yoshizawa & Noguchi (2003) 
included, for the first time, gas dynamics and star forma- 
tion processes, while representing the SMC also as a rota- 
tionally supported exponential disk with a nearly spher- 
ical halo, having the disk-to-halo mass ratio of 3:7. The 
significantly lower disk mass resulted in the disk being sta- 
ble against bar instabilities. Following typical findings for 
the Magellanic- type galaxies, a slowly rising rotation curve 
was assumed in both models, with a turnover radius of 3.5 
kpc and maximum rotation velocity of 50 km s^^. The 
best spatial orientation of the original SMC disk was de- 
termined in Gardiner & Noguchi (1996) to have an inclina- 
tion i = 45° and major axis 9 = 230° in order to match the 
observed gas and stellar distributions. These values were 
adopted by Yoshizawa & Noguchi (2003) and simulations 
were repeated for different star formation parameters. 

The major and common result of all three models is that 
the current 3-D structure of the SMC is composed of a 
central, disk- or bar-like, component and two tidal, spiral- 
arm-like tails. These tidal tails extend into the Magellanic 
Bridge, were formed 200 Myr ago, and are seen in both 
gas and stellar components. In particular, simulations 
by Yoshizawa & Noguchi (2003) show great morpholog- 
ical change of the initial SMC's gas disk - from its original 
size of about 10 kpc in diameter the disk shrunk forming 
first the Magellanic Stream and the Leading Arm (about 
1.5 Gyr ago), and then later two tidal arms that form the 
Magellanic Bridge (about 200 Myr ago) . At the end of the 
simulation the gas disk is still present but is now signif- 
icantly smaller in size, approximately by a factor of 2.5. 
The current SMC consists of the following components. 

1. The left-over disk- like component, located at 
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0*^15™ <RA< l'\ is not greatly elongated along 
the line-of-sight, with the stellar disk having a 
higher distance dispersion than the gas disk. 
Kinematically, the gas disk has a significant 
velocity gradient, from ^80 to 180 km s"-'^. The 
left-over stellar disk, on the other hand, shows high 
dispersion in the velocity field but no significant 
velocity gradient. In simulations by Gardiner & 
Noguchi (1996), this is a bar-like feature slightly 
bigger and centered more westward, covering 
O'^aO™ <RA< l^^aO™, and a slightly wider velocity 
range. 

2. The eastern tail, starting around RA ^ l'\ 
extending into the Magellanic Bridge, and covering 
a distance range from 55 to 40 kpc. This feature 
starts at a heliocentric velocity of about 170 

km s"""^ and gradually decreases to about 120 
km s~^ at RA > 2^. It is seen in both gas and 
stars but the gas tail is more extended. 

3. The western tail, starting around RA ~ 0.5^ and 

covering a distance range from 55 to almost 80 
kpc. Kinematically, this feature starts westward 
from the main disk-like component at a velocity of 
about 80 km s~^, but then turns east passing the 
main component with increasing velocity all the 
way to about > 250 km s"^ around RA > 2^. The 
stellar western tail follows the gas one and is less 
extended. 

Hence, the two tidal tails contribute mostly to the SMC's 
elongation along the line-of-sight, while the SMC disk- like 
component is not significantly elongated ( ~5 kpc). 

The best simulation by Yoshizawa & Noguchi (2003) re- 
produces a large number of observational properties in the 
Magellanic System, suggesting that most of them may be 
predominantly of tidal origin. Very importantly, a good 
morphological and kinematic reproduction of the Magel- 
lanic Stream with almost no stars was achieved for the first 
time, by using a very compact initial configuration of the 
SMC's stellar disk. The major shortcoming of the model 
is its failure to reproduce gas masses and mass ratios in 
the Magellanic Stream, Bridge, and the SMC. This may 
be related with the need for a more massive initial disk 
of the SMC, or a different halo-to-disk mass ratio, and re- 
quires further investigation. In addition, as pointed out 
by Putman et al. (2003), a number of detailed features 
related to the Magellanic Stream (e.g. its double-helix ap- 
pearance and numerous small clouds surrounding its main 
filaments) need to be explained. It has also been shown 
that forces other than gravity must play a significant role 
in the kinematics of the SMC (Kunkel et al. 2000; Zaritsky 
et al. 2000). 

The results of tidal numerical simulations show that the 
current central gaseous body of the SMC (^ 4 degrees 
in extent) has a significant velocity gradient. Gardiner & 
Noguchi (1996) explained this gradient as being due to the 
elongation of the SMC's bar-like component along the line- 
of-sight. However, a significant angular momentum left 
from before the last two close encounters with the LMC 
is expected to be present, if the pre-encounter SMC disk 
had an internal spin. The N-body simulations by Mayer 



et al. (2001) showed that in general disk-like dwarf irreg- 
ular galaxies undergo a great morphological transforma- 
tion due to close encounters with the more massive Milky 
Way. While dwarfs lose their mass through tidal stripping, 
their angular momentum gets gradually removed over a 
period of 6-10 Gyrs. This transformation is gradual, and 
on short timcscales (~ 2 Gyrs) a significant angular mo- 
mentum (>80%) is still preserved. Focusing especially on 
the SMC similar conclusion was reached by Kunkel 

et al. (2000) who found that it is extremely hard for the 
SMC to lose its pre-encounter angular momentum during 
an encounter with the LMC. 

5.2. Comparison with HI 

In Section 3.2 we showed that the HI distribution con- 
tains a large velocity gradient and some signatures for the 
existence of ordered, systematic motions. This agrees with 
theoretical simulations which show that the current cen- 
tral gaseous body of the SMC should still posses angular 
momentum left over from before the last two close encoun- 
ters with the LMC. This angular momentum should be less 
significant in the case of older stellar populations, as in- 
deed has been found from observations (see Section 3.3 for 
a summary) . We now compare RA- velocity slices through 
the HI data cube, shown in Fig. 6, with the theoretical 
predictions, shown in Fig. 17 of Yoshizawa & Noguchi 
(2003) and Fig. 10 of Gardiner & Noguchi (1996). The 
area enclosed in our observations contains only the central 
region of the SMC (4.5° x 4.5°). 

Images in Fig. 6 show coherent large-scale features with 
a wealth of small-scale structure. Many of the features 
have already been interpreted as expanding shells of gas 
(Staveley-Smith et al. 1997; Stanimirovic et al. 1999). 
At a low Dec of ~ —74°, the two most dominant features 
are supergiant shells 37A (RA Ol'^ 25™) and 494A (OO'' 
35™). These shells were investigated in Stanimirovic et 

al. (1999). Around Dec 73°20', Fig. 6 shows an HI 

emission, bar-like, feature stretching over the RA range 
QQh 4Qm r-^oO*^ 55™, and having a velocity gradient from 
'^120 to ^^170 km s~^. This feature continues eastward 
as a complex network of filaments and shell-like features, 
culminating around Dec ~ — 72°37' with two long fila- 
ments centered at heliocentric velocities of 130 and 170 
km s~^. The coherent appearance of these features traced 
throughout the data cube was interpreted as the low-Dec 
hemisphere of the supergiant shell 304A (sec Stanimirovic 
et al. 1999). At a very high Dec of 71°55', a high- 
velocity component is apparent, at about 200 km s~^. 

We do not find obvious features that could correspond 
to the predicted eastern and western tails within the re- 
gion covered by our observations. To search for potential 
tidal tails running out of the central SMC main body an 
investigation of a larger area is essential; this will be pos- 
sible in the near future with new observations by Briins 
et al. (2002) and MuUcr et al. (2003), but is beyond the 
scope of this paper. Concerning the position and veloc- 
ity range of the central component, our observations are 
in better agreement with the Gardiner & Noguchi (1996) 
model predictions: the coherent structure seen from RA 
00*^ 30™ to 01*^ 30™ corresponds to the area where the 
central bar-like feature predicted by the model should be 
found, and its velocity span from 100 to 200 km s~^ is 



8 



S. Stanimirovic et al. 



Fig. 6.— [fe.jpeg] Ten RA-velocity slices through the combined ATCA and Parkes HI data cube of the SMC, starting with Dec -74°42'12" 
and ending with Dec -72°17'42" (J2000). The grey-scale range is: 0-45, 0-45, 0-55, 0-75, 0-75, 0-75, 0-75, 0-55, 0-45, 0-45 K respectively, 
with a linear transfer function. 



Fig. 7. — [f7.jpeg] continued from Fig. 6 



very close to the predicted velocity range for the bar-like 
component. The Yoshizawa & Noguchi (2003) model pre- 
dictions for the left-over disk component cover a smaller 
area and velocity range than wliat is found in the obser- 
vations. 

6. ROTATION CURVE AND MASS MODELING OF THE SMC 

As discussed in the previous section, the central 4° x 4° 
of the SMC, covered by the HI observations, most likely 
corresponds to the left-over part of the original SMC's 
gaseous rotationally supported disk. This disk should con- 
tain significant signatures of angular momentum left from 
the pre-encounter disk. Motivated by this, we proceed 
further to derive the rotation curve from the observed HI 
velocity field, and the mass model, assuming that the en- 
tire velocity gradient is due to rotation. 

6.1. Tilted ring analysis 

The velocity field derived in Section 3.2, and shown in 
Fig. 5, was used to define the major and minor kinematic 
axes, and the apparent kinematic center at ^ RA Ol'^OS™, 
Dec -72°25' (J2000), with a systemic velocity of ~ 20 
km s~^ in the galactocentric reference frame. The po- 
sition angle (PA) of the major kinematic axis is around 
50°. The tilted ring algorithm ROCUR in the AiPS pack- 
age was then used to derive the HI rotation curve at dif- 
ferent radii from the center. A fairly standard procedure 
(Meurer et al. 1998; Bureau & Carignan 2002) was applied 
to estimate a set of kinematic parameters that represents 
well the observed velocity field at all radii. It was found 
that solutions for the inclination vary greatly especially 
for the outer radii. In the end, the value i = 40° ±20° was 
adopted. The newly determined least-squares solution for 
the systemic velocity was found to be 24 km s~^. Finally, 
the solutions for the rotational velocity and position angle 
were determined and are shown in Fig. 8. The position 
angle varies systematically from 30° to 50° (as shown in 
the top panel of Fig. 8). Using the mean PA, the global 
and the separate rotation curves for the receding and ap- 
proaching sides of the velocity field were derived, as shown 
in Fig. 8 (bottom panel). 

The global rotation curve shows a slow rise up to 
i? ~ 3 kpc, where it reaches the maximum velocity of ^^40 
km s~^. The receding and approaching curves are signifi- 
cantly different, which is not surprising since the velocity 
field is quite asymmetric (Fig. 5). The receding curve is 
less smooth and could have a local maximum around 1 
kpc from the center of rotation. The maximum rotation 
velocity found in previous studies (Hindman 1967; Loiseau 
& Bajaja 1981) was ~ 36 km s~^. 

The line-of-sight HI velocity dispersion distribution 
shown in Fig. 4 pointed to a high velocity dispersion across 
most of the SMC. The mean value, cthi = (22 ±2) km s"""^, 
is significantly higher than what is found for spiral galax- 
ies, or even other dwarf galaxies (Meurer et al. 1998; Cote 



et al. 2000). It is also large compared to the rotational 
velocity (Fig. 8). If interpreted as being all due to random 

motions, rather than bulk motions along the line-of-sight, 
such high values for cthi suggest that the turbulence in the 
ISM of the SMC has an important influence on the system 
dynamics. 

To account for this dynamical support we estimated the 

asymmetric drift correction. The usual prescription for 
determining this correction (Meurer et al. 1996) was fol- 
lowed: we derived the azimuthally averaged HI surface 
brightness and velocity dispersion profiles (for deprojected 
circular annuli with a mean position angle of 40° and an 
inclination of 40°), and assumed a constant vertical scale 
height. In general, the derived asymmetric drift correc- 
tion is quite significant and range from 0, for the inner 
radii, to ~ 40 km s~^ for the outer radii. The observed 
rotation curve corrected for the pressure support is pre- 
sented in Fig. 9 (crosses). The difference is significant for 
i? > 0.5 kpc with the corrected curve being higher by ~ 
10 km s~^ than the observed one. 

6.2. The mass model 

A two-component mass model was fitted to the corrected 
rotation curve. 

(a) To derive the deprojected rotation curve due to the 
potential resulting from the neutral gas alone, the radial 
HI surface density profile was used in the gipsy's rotmod 
task (van der Hulst et al. 1992). An exponential density 
law was assumed for the vertical disk distribution, with 
a scale height of 1 kpc. An indication of the large scale 
height was given by the typical size of large HI shells in the 
SMC, as well as from previous estimates based on the ve- 
locity dispersion and an average surface density of matter 
for the case of an isothermal disk. The resultant rotation 
curve {Vg) is shown in Fig. 9 with dotted line. 

(b) To derive the deprojected rotation curve due to 
the optical surface density distribution, the V-band im- 
age (shown in Fig. 2) was used. This image was first 
smoothed to 20 arcmin resolution and scaled to match the 
total luminosity of the SMC in the V-band of 3.1 x 10^ L© 
(de Vaucouleurs et al. 1991). The rotation curve arising 
from the stellar potential alone (14) was estimated from 
the V-band stellar density profile assuming, in first ap- 
proximation, that M*/Lv = 1 (see the curve shown as 
a dash-dotted line in Fig. 9). The best fit of the total 

predicted rotational velocity, -^V^ + V^, to the observed 

rotational velocity corrected for the pressure support was 
obtained for M^/Ly = 1.02 (solid line in Fig. 9). The two 
component mass model appears to fit the observed rota- 
tional velocities quite well, suggesting that no additional 
component, such as a dark halo, is needed to explain the 
rotation of the SMC. 

The rotation velocity of only the stellar component im- 
plies a total stellar mass of the SMC of 1.8 x 10^ M©, 
within a radius of 3 kpc. Prom Section 2 and after correc- 
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Fig. 8. — The least-squares solution for the position angle (PA) of the kinematic line of nodes (top panel) and the rotational velocity (V^ot) 
(bottom panel) as a function of radius. The mean position angle is 40° . The rotation curves for the approaching and receding sides of the 
velocity field are shown with filled and open circles, respectively. The global rotation curve is shown with the solid line. 
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tion for neutral Ho, the mass of HI+He is 5.6 x 10* M©, 
being almost a tliird of the stellar mass within the same 
radius. The total mass of the SMC, implied from the ro- 
tation curve, is thus 2.4 x 10^ Mq. This is almost twice 
that of 1.5 X 10^ Mq, derived by Hindman (1967) within 
the slightly smaller radius of 2.6 kpc, and is similar to the 
inital mass of the SMC assumed by Gardiner et al. (1994) 
(about 2 X 10-* M0) and by Yoshizawa & Noguchi (2003) 

(3 X 10^ Mq). 

6.2.1. The stellar M/L from population synthesis models 

Stellar population synthesis (SPS) models provide an 
independent way to estimate the stellar mass-to-light ra- 
tio. Bell & de Jong (2001) used simplified spiral galaxy 
evolution models, based on several different SPS models, 
to investigate trends of the stellar M/L with galaxy prop- 
erties, such as colors and gas fraction. They found good 
agreement with observational data and showed that the 
trend of M/L with color, logiQ{M/L) = a\ + ^Adolor, is 
robust with respect to the SPS models but is dependent 
on the type of the initial mass function (IMF) used. They 
tabulate coefficients a a and a a for different colors, stellar 
metalUcity, SPS models, and IMFs. 

We have used coefficients for the case of a stellar mctal- 
licity appropriate for the SMC, which is about 25% of the 
solar abundance, based on the SPS model by Bruzual & 
Chariot (2003), for a Salpeter IMF scaled by a factor of 
0.7 (Bell & de Jong 2001), and the Schmidt (1959) star 
formation law. The stellar M/Ly was then estimated us- 
ing the broad band colors given in Table 1, B — V = 0.41 
and B-R= 0.7. We get M/Ly = 0.8 ± 0.2 which agrees 
with the value obtained from the HI rotation curve. 

6.3. Critical density for star formation 

We now go a step further and address the gravitational 
stability of the current gaseous disk of the SMC and its 
possible consequences for star formation. We estimated 
the disk stability parameter, Q, (Toomre 1964) using the 
mean velocity dispersion, azimuthally averaged HI surface 
brightness, and the angular velocity derived using the ro- 
tation velocity. According to the Toomre's stability crite- 
rion, the disk is stable for Q > 1, while regions with Q ^ 1 
are unstable and under possible fragmentation which can 
lead to new star formation (Binney & Tremaine 1987). 
The estimated values of Q, as a function of radius, are 
shown in Fig. 10. The plot shows that most of the gaseous 
disk, for R < 2.1 kpc, is unstable and, most likely, under 
current star formation. The unstable region appears to 
be quite extended, reaching almost the cut-off Ha radius, 
i?(Hc) « 2.6 kpc (based on Ha observations by Kennicutt 
et al. (1995)), which delineates the parts of the gaseous 
disk under the most recent star formation. 

We also determined the so called critical surface den- 
sity, defined as the surface density needed to stabilize 
the disk (or for Q = 1). Both the observed (Shi) and 
critical (Sc) surface densities are shown in Fig. 11 (top 
panel). The ratio of the observed to critical surface density 
(a = Shi/Sc = 1/Q) is shown in the same figure (bottom 
panel) . For the unstable part of the disk, a varies between 
1.0 and 2.0, while it has a constant value of a = 0.74±0.04 
for the stable part of the disk. This agrees with the empir- 
ical relationship between the star formation threshold and 



the disk stability found by Kennicutt (1989) for a sample 
of normal galaxies, whereby most of the star formation is 
taking place for a > 0.67 (and hence for Q < 1.5). 

7. SUMMARY AND CONCLUDING REMARKS 

We have used the latest HI observations, obtained with 
the ATCA and the Parkes telescope, to re-examine the 
kinematics of the SMC. The HI velocity field, derived in 
Section 3.2, shows a large velocity gradient from the south- 
west to the north-east. The iso-velocity contours of this 
velocity field show some symmetry, suggestive of a differ- 
ential rotation. Some large-scale distortions in this veloc- 
ity field are easily visible, but could be related to posi- 
tions of several supergiant shells. In contrast to the HI 
distribution, the old stellar populations appear to have a 
spheroidal spatial distribution and a total absence of rota- 
tion. In Section 4, we summarized previous observational 
models of the SMC and the greatly controversial issue of its 
possibly large depth along the line-of-sight. We cautioned 
that the differential extinction across the SMC, revealed 
recently by Zaritsky et al. (2002), could be a significant 
source, unappreciated previously, for an overestimate of 
stellar distances, and therefore the SMC's depth. 

In Section 5, we summarized several tidal models that 
are concerned with the 3-D structure and kinematics of the 
SMC. In order to reproduce observational characteristics 
of the Magellanic System these models needed to assume 
the existence of angular momentum in the pre-encounter 
SMC (about 1.5-2 Gyr ago). Models predict that the cur- 
rent, left-over material from the pre-encounter SMC disk is 
a disk-like or a bar-like feature, about 4 kpc in extent, with 
a velocity gradient of about 100 km s^^. This velocity gra- 
dient is at least partially from the original angular momen- 
tum, as it is very difficult to lose the original spin during a 
galaxy encounter (Mayer et al. 2001; Kunkel et al. 2000). 
The elongation along the line-of-sight of the pre-encounter 
SMC disk may also be partially responsible. It is, however, 
not clear from the models what are the relative contribu- 
tions of these two effects to the predicted velocity gradient. 
In addition, theoretical models predict higher velocity dis- 
persion for the post-encounter stellar disk without a sig- 
nificant velocity gradient; this agrees with observations of 
older stellar populations summarized in Section 3.3.2. 

While we do not find evidence for the existence of tidal 
tails predicted by the models, the central region of the 
SMC covered in these observations most likely corresponds 
to the central component left-over from the pre-encounter 
SMC gaseous (rotationally supported) disk. The observed 
HI velocity gradient agrees well with the model predic- 
tions. We then proceeded further to derive its rotation 
curve and mass model in Section 6. A set of kinematic 
parameters derived from the tilted ring analyses agrees 
extremely well with the assumptions used in theoretical 
models that led to a good reproduction of observational 
properties of the Magellanic System. The derived HI rota- 
tion curve in Section 6 rapidly rises to about 60 km s~^ up 
to the turnover radius of ~ 3 kpc. A stellar mass-to-light 
ratio of about unity was required to scale the stellar com- 
ponent of the SMC's rotation curve to match the observed 
rotation curve. This suggests that a dark matter halo is 
not needed to explain the dynamics of the SMC. This is 
a surprising result, as dwarf irregular galaxies are often 
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Fig. 11. — Top panel: the radial distribution of the observed (Shi) and critical (Sc) surface density of the gas in the SMC. Bottom panel: 
the radial distribution of the ratio of the observed to the critical surface density (a = Ehi/Sc). The dashed line shows a = 0.74. 
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found to be dark matter dominated, and could be related 
to the mechanism of tidal stripping, as indicated by Mayer 
et al. (2001). The total dynamical mass of the SMC de- 
rived from the rotation curve is 2.4 x 10^ M©, three quar- 
ters of which are due to the stellar mass only. Wc also 
derived Toomre's disk stability parameter Q, which shows 
that almost all of the SMC disk is in the unstable regime. 
The gravitationally stable part, on the other hand, has a 
ratio of the observed to critical surface density a = 0.74 
which is in excxiknit agreement with the empirical stability 
threshold found by Kennicutt (1989). 

All of the above suggest that the HI distribution is ca- 
pable of providing valuable information about the SMC 
dynamics. Although very disturbed, it still contains im- 
prints of the original system prior to the latest encounters 
with the LMC and the Galaxy. The HI data discussed 
here comprise only the central 4.5° x 4.5° of the SMC. It 
is very important to compare tidal theoretical models with 
data comprising a larger area around the SMC to search 



for signatures of the proposed tidal tails. This will be pos- 
sible in the near future using the Parkcs Multibeam HI 
observations of the whole Magellanic System (Briins et al. 
2000), as well as the ATCA and Parkes observations of the 
Magellanic Bridge (MuUcr et al. 2003). 

The HI SMC data set is available from the ATNF SMC 
web page 

(http: / /www. atnf . csiro.au /research /smcJil /) . 
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